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Abstract 


Changes in precipitation patterns and its intensity due to the global climate change has led to 
intensifying floods as the most occurent natural disaster. Flooding cannot be fully prevented, but 
its impacts can be mitigated by accurately identifying flood-prone areas and implmenting 
appropriate risk-based management measures. The Qarachai Watershed, located in the upstream 
of the Gorganrood River Basin in Golestan Province, Iran, was chosen as the study area, which 
has experienced several flood events in recent decades. In this study, flood risk was assessed using 
the HEC-RAS, a two-dimensional hydraulic model. Flood discharge values occured in the early 
2019 were considered as model input, and based on ground truth, the manning roughness 
coefficient values were measured. To evaluate the results of the HEC-RAS model, some statistical 
criteria were used reflecting a good performance of the model. The analysis showed that by 
increasing the return period, the extent, depth, and amount of flood risk increase. Moreover, the 
analysis showed that the flood zone in 100 year return period, affects parts of Seyedkalate Village. 
Approximately, half of the flood zone identified in this study was attributed with very low-risk. 
The results of the study are used to adopt appropriate strategies and plans to adapt to climate 
change and as an appropriate tool for identifying flood exposed and flood-prone zones. 
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حغرافیا و مخاطرات محیطی. سال دوازدهم. شمارة جهل و هشتم. زمستان ۰۱۶۰۲ صص ۱۸۷-۲۰۳ 
مقاله پژوهشی 


بررسی قابلیت مدل دویعدی HEC-RAS‏ در تهیه نفشه خطر سیلاب در رودخانه قره‌جای رامیان. 
استان گلستان 


محبوبه بای- دانشجوی دکتری مهندسی آبخیزداری دانشگاه لرستان, = able‏ ایران 
8 ناصر طهماسبی پو !— دانشیان هیدرولوژی, دانشگاه oly! ble = Oki JJ‏ 
حسین زینی‌وند- دانشیار» هیدرولوژی. دانشگاه لرستان» = ble‏ ایران 
امیر سعدالدین = استاد. مدیریت آبخیز دانشگاه علوم‌کشاورزی و منابع‌طبیعی گرگان. گرگان ایران 
مهدی کاهه- دکتری منابع آب (سازه‌های هیدرولیکی» شرکت سهامی آب منطقه‌ای گلستان» گرگان» ایران 


تاریخ دریافت: ۱2۰۰/۱۲/۸ تاریخ بازنگری: ۱8۰۱/۲/۲۵ تاریخ تصویب: ۱۶۰۱/۳/۶ 


حکیده 


تغییرات در الگوها و شدت‌های بارش ناشی از تغییر اقلیم جهانی موجب تشدید رخدادهای سیل به‌عنوان 
یکی از رایج‌ترین بلای طبیعی در سطح دنیا شده است. این در حالی است که بشر نمی‌تواند جلوی بروز 
سیل را بگیرد. اما می‌تواند اثرات این مخاطره را با شناسایی دقیق مناطق مستعد سیل و اجرای اقدامات 
مدیریت ریسک محور کاهش دهد. حوزه آبخیز قره‌چای یکی از سرشاخه‌های حوضه گرگانرود در استان 
گلستان به‌عنوان منطقه موردمطالعه انتخاب شد که رخدادهای سیلابی متعددی را تجربه نموده است. در 
این مطالعه خطر سیل با استفاده از مدل هیدرولیکی HEC-RAS‏ دوبعدی ارزیابی شده است. بدین‌منظور 
مقادیر دبی رخداد سیلابی سال ۱۳۹۷ به‌عنوان ورودی مدل مذکور انتخاب شدند و بهینه‌سازی ضرایب 
زبری مانینگ مطابق با شواهد موجود انجام شد و سپس واسنجی مدل انجام شد. شاخص‌های Solel‏ مورد 
استفاده در ارزیابی مدل همگی توصیف‌کننده شبیه‌سازی خوب آن می‌باشند. نتایج نشان داد با افزايش دوره 


باز گشست سیل. گستره عمق و میزان خطر ناشی از OT‏ افزایش می‌یابد؛ به‌طوری‌که با در نظر گرفتن دوره 


Email: Tahmasebi.n @lu.ac.ir ۰۹۱۲۵۷۸۷۳ نویسنده مسئول:‎ ۱ 

نحوه ارجاع به این مقاله: 

HEC- محبوبه؛ طهماسبی پور» ناصر؛ زینی وند. حسین؛ سعدالدین, امیر؛ کاهه» مهدی. (۱۶۰۲). بررسی قابلیت مدل دو بعدی‎ ch 

.۱۸۷-.۰۲۰۳)۶(۱۲ قره‌چای رامیان» استان گلستان. جغرافیا و محاطرات محیطی.‎ leas) نقشه خطر سیلاب در‎ agi در‎ RAS 
https://doi.org/10.22067/geoeh.2022.75557.1185 


3A۸‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


بازگشت صدسال در مطالعه حاضر. پهنه سیلاب بخش‌هایی از روستای سیدکلاته را تحت‌تأثیر قرار داده 
است. در نقشه پهنه‌بندی خطر سیل. حدود نیمی از بازه موردمطالعه در کلاس بسیار کم‌خطر (کلاس یک) 
قرار دارند. نتایج مطالعه به‌منظور اتخاذ راهبردها و برنامه‌های مناسب با هدف سازگاری با شرایط تغییر 
اقلیم و به‌عنوان ابزار مناسب برای شناسایی مناطق مستعد و در معرض خطر سیل کاربرد دارد. 

کلیدواژه‌ها: نقشه خطر سیل. مدل HEC-RAS‏ دوبعدی. مدیریت ریسک سیل, رودخانه قره‌چای, استان 
گلستان. 


۱- مقدمه 

در سال‌های اخیر تحقیقات اقلیمی گسترده‌ای به‌منظور پیش‌بینی تغییرات الگو و شدت بارش Ske‏ تا پایان قرن 
۱ انجام شده است. نتایج مطالعات حاکی از افزایش بارندگی در مناطق مختلف دنیا اسست. تغییرات ناگهانی در 
الگوهای بارش و شدت آن منجر به بروز خطرات طبیعی از قبیل سیل. خحشکسالی» زمین‌لغزش ناشی از بارش 
می‌شود. ازاین‌بین سبلاب‌ها بیشترین فراوانی و بالاترین خطر" را به خود احتصاص داده است (رنجری و همکاران ۰ 
۵۹ و مطالعات انجام شده در زمینه تغییر اقلیم بر افزایش تعداد افراد در معرض سیل دلالت دارند (تلمن و 
همکاران ؛ ۲۱ به‌طوری‌که در طول ۲۰ سال گذشته. سیل رایج‌ترین بلای طبیعی محسوب می‌شده است و EY‏ 
درصد از وقوع کلیه مخاطرات طبیعی (جوهااسپر و همکاران؛ ۲۰۱۵؛ دی‌بلداساری و همکاران؟ ۲۰۱۰) مربوط به آن 
است. سیل‌ها یکی از عوامل بالقوه مرگ و میر مردم و جابجایی آن‌ها از مساکن خود در سطح جهان هستند که ضمن 
آسیب به محیط‌زیست. اقتصاد توسعه اقتصادی آن منطقه را pees‏ قرار می‌دهد (دی‌بلداساری و همکاران. 
۰ درحالیکه نمی‌توان جلوی وقوع سیل را گرفت. اما اثرات این مخاطره را می‌توان با شناسایی دقیق مناطق 
مستعد سیل و اقدامات مدیریتی ریسک محور کاهش داد. مدیریت ریسک محور سیل به‌عنوان روشی مناسب برای 
مواجه با مخاطرات محیطیء محیطزیستی و خطرات ایجاد شده توسط انسان تعریف می‌شود CEL)‏ ۲۰۰۲). شناخت 
فرایندهای سیلابی شدن بر حسب اندازه سیل. عمق طغیان» سرعت جریان و زمان رسیدن سیل در تدوین 
استراتژی‌های کاهش خطر سیل ضرورت دارد (ساهو و سجا“ ۲۰۱۷)؛ بنابراین اتخاذ هر نوع اقدامی نیاز به کسب 


1 Hazard 

2 Rangari et al. 

3 Tellman et al. 

4 Guha-Sapir et al. 

5 Di Baldassarre et al. 
6 Plate 

7 Sahoo and Sreeja 

8 Amrei and Britta 


مدیریت موثر سیلاب نیازمند نقشه طغیان" سیل» تخمین احتمالی از آسیب و ریسک سیل در پهنه سیلاب و ارائه 
طرح جامع برای کاهش ریسک سیل است. نقشه خطر یک جزء مهم در مدیریت کارآمد ریسک به دلیل فراهم‌کردن 
اطلاعات عمق و پهنه سیل محسوب می‌شود (شهیری‌پارسا و همکاران » ۲۰۱۹). 

به‌منظور بررسی رفتار و ویژگی‌های سیل برای دشت‌های سیلابی و خطرات آن‌ها دو نوع رویکرد متمایز وجود 
دارد: (il‏ استفاده از سامانه اطلاعات جغرافیایی» تصاویر ماهوارهای و داده‌های تاریخی ثبت شده (داغاب و علائم 
سیل) که در مواقم کمبود داده و عدم وجود آن مورد استفاده قرار می‌گیرد. ب) استفاده از مدل‌های عددی برای 
Gilt‏ جریان‌های سیلابی که مورد اخیر به دلیل توانایی بیشتر در ارائه نتایج خوب. قابل قبول‌تر اسست و 
انعطاف‌پذیری آن باعث شده رفتارها و ویژگی‌های سیل به‌وضوح مشخص و بررسی شود (ترین و ملکنتین " ۲۰۲۱). 

از سویی توسعه مدل‌های هیدرولیکی همیشه یک ابزار و گام ضروری در آنالیز سیلاب‌ها می‌باشند. در سال‌های 
Gt!‏ عمدتاً از مدل‌های هیدرولیکی یک‌بعدی و دوبعدی به‌عنوان ابزاری برای ارزیابی سیل و ارزیابی خطر سیل 
استفاده می‌شود. اونقداز و همکاران؛ ۲۰۲۰ در رودخانه یسیل (اشیم)" قزاقستان اقدام به شبیه‌سازی پهنه سیلاب با 
مدل دوبعدی HEC-RAS‏ نمودند. نتایج مطالعه آن‌ها نشان داد تغییر اندازه شبکه‌ها در عملکرد مدل تغییر چندانی 
ایجاد نمی‌کند. همچنین بر اساس نقشه‌های تهیه شده پهنه‌بندی سیلاب با دوره بازگشت‌های مختلف موردمطالعه و 
طبقه‌بندی مناطق خطر سیل. مناطق مسکونی در بخش‌هایی از پایتخت قزاقستان در تمامی دوره بازگشت‌ها در زیر 
طغیان سیل قرار داشستند. امری و بریتا (۲۰۲۰) طی مطالعه‌ای در حوزه آبخیز فیش‌باخ" در جنوب آن به مقایسه 
رویکردهای هیدرولوژیکی و هیدرودینامیکی (دو رویکرد یکپارچه و جداشده) به‌منظور تهیه نقشه خطر سیلاب 
پرداختند. هیدرولوژی حوزه آبخیز با استفاده از مدل‌های نیمه‌توزیعی بارش رواناب Hec HMS‏ تعیین و نقشه خطر 
سیل با استفاده از هیدروگراف‌های زیرحوزه‌ها در یک مدل هیدرودینامیکی HEC- RAS‏ دوبعدی تهیه شد. نتایج 
مطالعه OLE‏ می‌دهد رویکرد جدا شده نتایج بهتری در هیدروگراف نسبت به رویکرد یکپارچه نشان می‌دهد. پوینز و 
تیم" (۲۰۱۹) در مطالعه‌ای در بازه رودخانه به‌اندازه پنج کیلومتر عملکرد چهار مدل هیدرولیکی بر اساس گستره و 
ارتفاع سیل به‌منظور شبیه‌سازی رخداد سیل با هدف تولید نقشه‌های مناطق سیل‌زده در رودخانه سانتا باربا" واقع در 


1 Inundation mapping 

2 Shahiri Parsa et al. 

3 Trinh and Molkenthin 

4 Ongdas et al. 

5 River Yesil (Ishim) 

6 Fischbach 

7 Decoupled and Integrated approach 
8 Pinos and Timbe 

9 Santa Barbara 


1۹۰ جغرافیا و مخاطرات محیطی شمارة چهارم 


آنداکوادور' در نزدیکی شهر گوالاستو" را مورد مقایسه قرار دادند. نتایج آن‌ها نشان داد در دوره بازگشت‌های ۲ و ۵۰ 
سال به ترتیب HEC-RAS 2D Flood Modeller 2D [ber 2D (clade‏ و PCSWMM 2D‏ عملکرد بهتری را 
نشان می‌دهند. رنجری و همکاران (۲۰۱۹) به ارزیابی ریسک آبگرفتگی در سیلاب‌های شهری با استفاده از 
مدل HEC- RAS‏ دوبعدی در حیدرآباد هندوستان پرداختند. نتایج نشان داد ۱۷ درصد از مناطق هدف موردمطالعه در 
معرض سیل قرار دارند که عمدتاً در کلاس ریسک متوسط قرار دارند. کومار و همکاران" (۲۰۱۹) از مدل HEC-‏ 
5 دوبعدی و سیستم جهانی پایش سیل " گستره سیل و نواحی با ریسک‌های مختلف در رودخانه گانا و "GUL‏ 
شهر پارایاگراجی" هند را بررسی کردند. نتایج OLS‏ داد که حروجی مدل HEC- RAS‏ دوبعدی و سیستم پایش 
Ske‏ با داده‌های مشاهداتی مطابقت دارد. آرنل و گاسلینگ" (۲۰۱۳) به بررسی اثرات تغییر اقلیم بر ریسک سیل‌های 
رودخانه‌ای در مقیاس جهانی با استفاده از چهار شاخص تغییر بزرگی و دوره بازگشت دبی اوج سیل. جمعیت 
موجود در مناطق سیل‌زده و زمین‌های زراعی در معرض سیل پرداختند. نتایج آن‌ها OLS‏ می‌دهد. تحت مدل تغییر 
اقلیم و سناریوهای انتشار گازهای گلخانه‌ای و اقتصادی اجتماعی در سال ۲۰۵۰ سیل‌های کنونی با دوره بازگشت 
۰ سال در 4۰ درصد سطح boo‏ حداقل دو برابر بیشتر اتفاق خواهد افتاد و40۰ میلیون نفر در مناطق سیل‌زده قرار 
می‌گیرند و 6۳۰ هزار کیلومترمربع از مزارع جزو مناطق سیل‌زده و در معرض دو برابر شدن فراوانی سیل خواهند بود. 
قنبرپور و همکاران" (۲۰۱۱) در مطالعه‌ای به واسنجی مدل هیدرولیکی HEC-RAS‏ براساس داده‌های مشاهداتی 
اندازه‌گیری شده برای رخداد سیل ناگهانی سال ۱۹۹۹ با دوره بازگشت ۲۰۰ سال و با استفاده از ضریب زبری 
مانینگ" پرداختند. نتایج مقایسه بین عمق آب رخداد سیل شبیه‌سازی شده و مقادیر بهینه ضریب زبری مانینگ 
همبستگی معنی‌داری را نشان می‌دهد. 

مرور منابع انجام شده همگی بر عملکرد مطلوب مدل HEC- RAS‏ دوبعدی به‌منظور ارزیابی ریسک سیل و نیز 
تهیه نقشه خطر سیلاب دلالت دارد. لذا در مطالعه حاضر نیز از مدل مذکور به‌منظور شبیه‌سازی سیلاب استفاده شد. با 
توجه با متأثر بودن استان گلستان همچون سایر نقاط کشور و جهان از رخدادهای سیلابی ناشی از تغییر اقلیم و 
تشدید خسارات ناشی از آن ارزیابی ریسک سیل و تهیه نقشه خطر در آن ضرورت دارد. ازاین‌بین حوزه آبخیز 


قره‌چای شهرستان رامیان به دلیل تجربه رخدادهای سیلابی رودخانه‌ای و ناگهانی متعددی به‌عنوان حوضه موردمطالعه 


1 Ecuadorian Andes 

2 Gualaceo 

3 Kumar et al. 

4 Global Flood Monitoring System (GFMS) 
5 Ganga and Yamuna 

6 Prayagraj 

7 Arnell and Gosling 

8 Ghanbarpour et al. 

9 Manning's roughness coefficient 


Red TE‏ بو iS‏ تاه Jet‏ سا ۱۳۹۳ در اسان اسان کرابم و CS‏ شاه ا 
ag‏ نقشه خطر و مناطق طغیان سیل یک ابزار مدیریتی مناسب برای شناسایی مناطق سیل‌زده و در معرض سیل 
محسوب می‌شود تا به‌منظور سازگاری با شرایط تغییر اقلیم و بروز رخدادهای سیلابی ناشی از آن» بتوان برنامه‌ها و 
اسع اا هان اسب انشا کر طا مهار اکان غ یا و kee‏ عا ها از doa‏ 
هیدرولیکی HECRAS‏ دوبعدی Gly‏ بازه‌ای از رودخانه قره‌چای رامیان استان گلستان پرداعته است. پرخلاف 
کاربردهای رایج برای انتشار سیل در بازه مذکون نسخه HEC-RAS‏ ۵.۰۷ چندان در ادبیات مورد استفاده قرار نگرفته 


است و قابلیت و پتانسیل مدل در شبیه‌سازی سیلاب مورد تجزیه‌وتحلیل قرار می‌گیرد. 


۲- مواد و روش 

۱-۲- موقعیت منطقه موردمطالعه 

حوزه آبخیز قره‌چای به مختصات ۲۳ ۳۸ ۰۵7 ۰۳ ۳۷ طول شرقی و ۰۲۲۶ ۵۵ NV EVE‏ ۵۵ عرض 
شمالی» یکی از سرشاخه‌های حوزه رودخانه‌ای گرگانرود در استان گلستان به مساحت ۲٤۸۱۷۵‏ هکتار به شمار 
می‌آید. حداقل و حداکثر ارتفاع آن به ترتیب ۲۸۵۰ متر و ۱۵۰ متر از سطح دریا متغیر می‌باشد (اداره کل منابع طبیعی 
و آبخیزداری استان گلستان ۱۳۸۲). میانگین بارندگی سالانه حوضه قره‌چای رامیان ۸۱۰/۳ میلی‌متر است. با توجه به 
تأثیرپذیری حوزه آبخیز قره‌چای از عوامل طبیعی و نیز عوامل انسانی از قبیل شیب زیاد منطقه. غالب بودن سازند 
حساس DES‏ تجربه وقوع بارش‌های نقطه‌ای c‏ عدم رعایت اصول صحیح کشاورزی در اراضی شیبدار 
کم‌بازده» بهره‌برداری نادرست از معادن» ایجاد خطوط ارتباطی بدون در نظرگرفتن اصول فنی در جاده‌سازی» استقرار 
سکونت گاه‌ها در مخروطه افکنه‌ها و حروجی حوزه‌های آبخیز, رخدادهای سیلاب رودخانه‌ای و ناگهانی به همراه 
خسارات قابل‌ملاحظه به وقوع پیوسته است (شکل ۱). 


1 Flash flood 


۱۹۲ جغرافیا و مخاطرات محیطی شمارة چهارم 
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۱-۲-۲- تهیه هیدروگراف سیل صدساله 

مقادیر بارش با دوره بازگشت صدساله از منحنی شدت. مدت و فراوانی ' به روش وزیری (وزیری»۱۳۹۳) 
استخراج و سپس از مدل "11150-11715 برای شبیه‌سازی فرایندهای هیدرولوژیکی جریان و تهیه هیدروگراف سیل با 
دوره بازگشت صدسال برای زیرحوزه‌های منتهی به بازه موردمطالعه براساس روش سرویس حفاظت منابع طبیعی " 
آمریکا و محاسبه شماره منحنی * استفاده شد. بر اساس شرایط حاکم در منطقه حوزه آبخیز قره‌چای به ۱۱ زیرحوزه 

۲-۲-۲- روندیابی حریان سیلابی با استفاده از مدل دوبعدی HEC-RAS‏ 

به‌منظور تعیین ژئومتری و ویژگی‌های هیدرولیکی سطوح دوبعدی جریان از لایه هندسی زمین" با استفاده از 
تصاویر ماهواره‌ای Woldview-¥‏ سال ۲۰۱۹ با قدرت تفکیک مکانی نیم‌متر برای بازه موردمطالعه به طول هفت 


1 Intensity Duration Frequency Curve (IDF) 

2 Hydrologic Engineering Center — Hydrologic Modelling System (HEC-HMS) model 
3 Natural Resources Conservation Service (NRCS) 

4 Curve Number (CN) 

5 Terrain layer 


کیلومتر نقشه DEM!‏ تهیه شد و در ادامه gY‏ مذکور با استفاده از تعدادی نقاط کنترلی و بازه‌های موردبررسی با 
دوربین نقشه‌برداری میدانی. صحت‌سنجی نقشه مذکور انجام شد. مقادیر دبی ساعتی ثبت شده رخداد سیلابی اسفند 
۷ و مرداد ۱۳۹۶ در ایستگاه هیدرومتری حوزه آبخیز قره‌چای به‌عنوان واقعه مشاهداتی به ترتیب برای مرحله 
واسنجی و اعتبارسنجی مدل" در نظر گرفته شد و سپس از مقادیر هیدروگراف agi‏ شده از حروجی مدل HEC-‏ 
5 برای دوره بازگشت صدسال مدل برای شبیه‌سازی جریان در دوره بازگشت مذکور استفاده شد. دوره 
بازگشت صدسال به‌عنوان دوره بازگشت استاندارد و پایه در محاسبات سیل و مدیریت ریسک سیل توسط بسیاری 
از مؤسسات و بسیاری از آژانس‌ها (اتحادیه مدیران امور سیلاب دولت" ۲۰۰۶) انتخاب و مبنای مطالعه حاضر قرار 
گرفت. 

۳-۲-۲- مدل دویعدی HEC-RAS‏ 

مدل دوبعدی 18٨-848‏ یکی از متداول‌ترین نرم‌افزارهای مدل‌سازی سیلاب در شبیه‌سازی هیدرودینامیکی 
است که به‌وسیله مرکز مهندسی هیدرولوژیکی ارتش آمریکا* تولید شده است. توسعه روش‌های عددی بیشتر و 
قلبل‌اعتمادتر قدرت محاسباتی کارآمدتر و تکنیک‌های جدید تهیه نقشه توپوگرافی منجر به افزایش استفاده از 
مدل‌های دوبعدی در شبیه‌سازی رخدادهای سیلابی شده است. به‌طوری‌که باوجوداینکه زمان لازم برای پردازش 
کمتر شده است. ولی از دقت مناسب برخوردار است. مدل‌های یک‌بعدی از کارایی و توانایی تولید فرایندهای درون 
کانال‌ها برخوردار هستند» Jy‏ مشکلات از زمانی شروع می‌شود که جریان سیلابی از کناره‌های JUS‏ سرریز* نموده 
و وارد دشت سیلابی می‌شود (پوینز و تیم ۶۹ و کاستیبل و همکاران"» ۲۰۲۰). از بین مدل‌های دوبعدی» مدل 
دوبعدی HEC-RAS‏ که از معادلات با مبنای فیزیکی " استفاده می‌نماید که یکی بهترین مدل‌های شناخته شده و 
پرکاربرد در ادبیات علمی و کاربردی محسوب می‌شود (کاستیبل و همکاران. ۲۰۲۰). همچنین به دلیل غالب بودن 
بعد افقی بر عمودی در Ol‏ مدل مذکور می‌توانند توصیف واقع‌بینان‌تری را از سیستم ارائه دهد (پوینز و تیم ۲۰۱۹). 
در مدل مذکور از معادلات دوبعدی سنت‌ونانت" یا معادلات موج پخشیدگی دوبعدی برای شبیه‌سازی جریان استفاده 
می‌کنند. به‌طورکلی معادلات موج پخشیدگی دوبعدی موجب اجرای مدل با سرعت بیشتری می‌شود و پایداری آن 


به‌مراتب بیشتر است. معادلات دوبعدی سنت‌ونانت دوبعدی برای محدوده وسیعی از مسائل قابل کاربرد هستند هر 


1 Digital Elevation Model (DEM) 

2 Model calibration and validation 

3 Association of State Floodplain Managers 

4 U. S. Army Engineer Hydrologic Engineering Center (HEC) 
5 Overbank flow 

6 Costabile et al. 

7 Physically- Based 

8 Saint-Venant's Sain 


۱۹۶ جغرافیا و مخاطرات محیطی شمارة چهارم 


چند در بسیاری از مدل‌سازی‌ها معادلات موج پخشیدگی در مقایسه با معادلات موج دیفیوژن نتایج دقیق‌تری را ارائه 
می‌دهند. در مدل دوبعدی HEC-RAS‏ این امکان وجود دارد تا بتوان با داده‌های در دس‌ترس و نوع پروژه موردنظر 
یکی از دو روش را انتخاب نمود (سیستم آنالیز رودخانه‌ای HEC-RAS .)۲۰۱۳ HEC-RAS'‏ یک ابزار پیچیده 
است که شامل تعداد زیادی برنامه برای مدل یک‌بعدی است و اخیراً معادلات آب‌های کم‌عمق " دوبعدی را به‌طور 


کامل که به شرح زیر نوشته شده است حل می‌کند Fey) ES‏ کاستیبل و همکاران Ye Ye‏ مویاکیوریگا و 
همکاران؛ ۲۰۱۲ و براثر ۲۰۱7). 


رابطه (۱) م 244+ 


0 a ۲2۶ , a 29 a fpepge aH ۳ ۳ 
۲+۶ (2) + 1 (24) = "Pave Loghi FPF + 55 (htxx) + 555 (htxy) )۲( رابطه‎ 


at 22۳۱۳ 
۹4 ۶ (pd) ب‎ 9 (a) = Eagt _ رن‎ 2 2 
+ )2( + 2 (£) = h2 9 برچ‎ t+ af + Fa, (xy) + رچم‎ (yy) )۳( رابطه‎ 


که در آن 6,۵ + Hyt) = zxy)‏ ارتفاع سطح Z (m)‏ ارتفاع سلول در مختصات کاتزین" × و Dg y‏ عمق 
آب <m)‏ ۷ط = و hu and‏ = و = UL‏ ویژه در جهت × و gu ۳02 8-1( y‏ ۷ به ترتیب سرعت‌ها در × و TY‏ شبکه 
باران ۲ am)‏ ع: شتاب تقل“ dms)‏ : ضریب زبری مانینگ dsm")‏ م چگالی dkg m°) “Ol‏ ××ے TXY y TYY‏ 
اجزای فشار تانسور" ؛ f‏ پارامتر کوریولیس (5-1). زمانی که موج دیفیوژن انتخاب شده است. 

شرایط مرزی بالادست و پایین‌دست در بازه موردمطالعه. به‌منظور تعیین مقادیر شرایط مرزی برای شبیه‌سازی 
جریان رودخانه غیرماندگار بر اساس رخدادهای مشاهداتی و هیدروگراف شبیه‌سازی از مدل HEC-HMS‏ برای 
دوره بازگشت تعیین شدند. از ضریب زبری مانینگ که یک پارامتر مهم هیدرولیکی در مدل 1150-8۸5 است که 
برای واسنجی مدل استفاده شد. در طول جریان» محاسباتی مربوط به سطوح دوبعدی مورد استفاده قرار گرفت تا 
مشخصات هیدرولیکی جریان در هر سلول به دست آید. به‌منظور ساخت نقشه توزیع مکانی ضریب زبری'' منطقه 
موردمطالعه از نقشه کاربری از تصاویر ماهواره‌ای لندست ۲۰۱۹ و ستتینل AY‏ ۲۰۱۹ در محیط GIS‏ استفاده 


1 HEC-RAS River Analysis System 
2 Shallow Water Equations 

3 Tabarak 

4 Moya Quiroga et al. 

5 Brunner 

6 Cartesian coordinate 

7 Net rain 

8 Gravity acceleration 

9 Water density 

10 Stress tensor 

11 Spatial Distribution of Manning Coefficient 


شد و مقادیر ضریب مانینگ کاربری‌های موجود طی بازدیدهای میدانی» بر اساس مطالعات go)‏ ۱۹۵۹ و آرسمنت 


و اشنایدر" ۱۹۸۹) و نیز تأیید کارشناسان و متحصصین تعیین شد (جدول ۱). مقادیر ضریب مانینگ برای رودخانه 


قره‌چای بر اساس تقسیم زیرحوزه به مقاطع طولی و به روش TOS‏ تعیین شد (رابطه AE‏ 


n=(no+ni+n2+n3+n4) M5 ٤ رابطه‎ 


که در آن 7 ضریب زبری مانینگ. 90 مقدار al‏ برای کانال صاف. یکنواختی و مستقیم در مواد طبیعی» nl‏ 
ضریب مربوط به آثار نامنظمی مقاطع عرضی, 2 چگونگی تغیبرات در اندازه و شکل MB SUS‏ وجود موانع در 
مسیر SUS‏ 94 نوع و تراکم پوششگیاهی» mS‏ درجه مثاندری بودن مسیر است (آرسمنت و اشنایدن .)۱۹۸٩‏ 

1-۲-۲ - واسنحی و اعتبارسنحی مدل 

در شبیه‌سازی جریان سیلابی ضریب زبری مانینگ رودخانه حساس‌ترین پارامتر محسوب می‌شود. برای 
شبیه‌سازی رودخانه با استفاده از مدل HEC- RAS‏ دوبعدی در بازه مورد نظر. مقادیر مذکور رودخانه و دشت‌های 
سیلابی لازم است با تنظیم مجموعه‌ای از پارامترها تابع هدف * بهینه‌سازی شوند (تبارک و همکاران» ۲۰۲۱). همچنین 
کمبود و عدم دسترسی به سوابق ثبت شده سیل» واسنجی مدل را با مشکل جدی مواجه می‌نماید. لذا بر اساس 
داده‌های داغاب» عمق سیل و سطح بالای تراز سیل انجام می‌شود (پوینز و تیم» ۲۰۱۹ و نعیم و همکاران؛ ۲۰۲۱). در 
مطالعه حاضر واسنجی مدل بر اساس بهینه‌سازی مقادیر ضریب زبری مانینگ براساس داغاب‌های موجود در منطقه و 
جمع‌آوری اطلاعات محلی مربوط به رخداد سیلاب ۱۳۹۷ انجام گرفت و با ثابت در نظر گرفتن ضرایب بهینه شده 
مانینگ, مدل بر اساس رخداد سال ۱۳۹۲ اعتبارسنجی شد. 

۵-۲-۲-رزیایی مدل 

به‌منظور ارزیابی مدل HEC-RAS‏ دوبعدی از شاخص‌های آماری مختلف استفاده می‌شود. شاخص‌های مورد 
استفاده در مطالعه حاضر عبارتند از: شاحص کارایی نش ساتکلیف": اهمیت نسبی واریانس مقادیر شبیه‌سازی شده را 
در مقایسه با واریانس داده‌های مشاهداتی OLE‏ می‌دهد و مقدار آن از منفی بی‌نهایت تا یک تغییر می‌نماید (رابطه ۵). 


اگر مقدار NSE‏ برابر با یک باشد نشان‌دهنده این است که تناسب کامل بین داده‌های مشاهداتی و شبیه‌سازی وجود 


1 Chow 

2 Arcement and Schneider 

3 Cowan 

4 Objective function 

5 Naeem et al. 

6 Nash-—Sutcliffe Efficiency (NSE) 


1۹1 جغرافیا و مخاطرات محیطی شمارة چهارم 


دارد» اگر مقادیر بین ۰/۳۳ تا ۰/۷۵ باشد. نتایج مدل رضایت‌بخش به شمار می‌آید و درصورتی‌که مقادیر بیش از ۰۷۵ 
باشد نتایج مدل خوب خواهد بود (مک‌اینتایر و القریشی! ۲۰۰۹). شاخص توافق" : شاخص مذکور توسط 
ویلموت (1۹۸۱) به‌عنوان معیار استاندارد شده برای تعیین میزان حطای پیش‌بینی شده مدل توسعه Cab‏ (رابطه (WV‏ که 
مقدار آن بین صفر تا یک متغیر است (موریسی " و همکاران ۲۰۰۷ و سولر و EOL Kae‏ ۲۰۰۷). هرچقدر مقادیر 
شاخص توافق به یک نزدیک‌تر باشند» عملکرد شبیه‌سازی مدل بهتر است (سولر و همکاران, ۲۰۰۷). ازآنجایی که هر 
دو شاحص کارایی نش‌ساتکلیف و شاخحص توافق بسیار تحت‌تأثیر مقادیر با ارزش‌های بالا هستند» به‌منظور کاهش 
خطا از معیارهای کارایی نسبی چون انحراف نسبی شاخحص کارایی نش ساتکلیف" (رابطه ۷) و انحراف نسبی 
شاخحص توافق" (رابطه (A‏ استفاده شد. انحراف نسبی شاحص کارایی نش‌ساتکلیف و انحراف نسبی شاحص توافق 
بین صفر و یک متغیر است. مقدار شاخص صفر و یک به ترتیب نشان‌دهنده عدم توافق و توافق است. (نعیم و 
همکاران. ۰۲۰۲۱ ریسو و همکاران" ۲۰۱۲؛ فقات و همکاران/ ۲۰۱7؛ کروس و همکاران ۲۰۰۵). 
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که در Qmean Ol‏ میانگین عمق آب مشاهداتی Om‏ عمق Ol‏ شبیه‌سازی شده QÉ‏ عمق آب مشاهداتی در 


1 McIntyre and Al-Qurashi 

2 Agreement 

3 Moriasi et al. 

4 Soler et al. 

5 Relative deviations of Nash—Sutcliffe efficiency (Erel) 
6 Relative deviations of Agreement (drel) 

7 Raposo et al. 

8 Phogat et al. 

9 Krause et al. 


۳- نتایج و بحث 


مقادیر ضریب مانینگ در طول جریان محاسباتی مربوط به سطوح دوبعدی برای تعیین مشخصات هیدرولیکی 
جریان در هر سلول مورد استفاده قرار می‌گیرد. جدول ۱ مقادیر ضریب زبری مانینگ کاربری‌های مختلف و بازه 
موردمطالعه رودخانه قره‌چای (هفت کیلومتر) نشان می‌دهد. همچنین نتایج مطالعه حاضر بر مبنای چهار شاخص 
otras E Slee J‏ شمان cl UY‏ بوسی ا داد انز تسب انا 
بین داده‌های مشاهداتی و شبیه‌سازی در مرحله واسنجی و نتایج رضایت‌بخش در مرحله اعتبارسنجی می‌شود و در 
بقیه شاخص‌های موردبررسی ald‏ بررسی شاخص‌های مورد استفاده همگی توصیف کننده شبیه‌سازی قابل‌قبول 
توسط مدل می‌باشد. 


جدول ۱- مقادیر ضریب مانینگ کاربری‌های مختلف در حوزه قره‌چای رامیان 


$ نام کاربری مقدار ضریب زبری مانینگ $ نام کاربری مقدار ضریب زبری مانینگ 
۱ اراضی زراعی NYY‏ ۶ بازه رودخانه۸ 0 
۲ دشت سیلابی a‏ ۷ بازه رودخانه٩ VPY‏ 
۳ جنگل f‏ ۸ بازه رودخانه ۱۰ ۰ 
۴ فعالیت کشاورزی AY‏ ۹ بازه رودخانه۱۱ PPA‏ 
۵ اراضی محدوده روستا Ys e‏ بازه رودخانه ۱۲ N‏ 
۶ مناطق مسکونی Neuss‏ ۳۱ بازه رودخانه۱۳ ANV‏ 
۷ مناطق اداری rade‏ ۳۲ بازه رودخانه۱۴ eA‏ 
۸ باغ vanes‏ ۲۳ بازه رودخانه۱۵ ELCT‏ 
4 بازه رودخانه۱ ۰ ۲۴ بازه رودخانه۱۶ ۰ 
۳ بازه رودخانه ۲ 1۰ Yò‏ بازه رودخانه۱۷ ۰ 
۱ بازه رودحانه۳ Af‏ ۲۶ بازه رودخانه۱۸ Vp‏ 
۱۲ بازه رودخانه۴ ANV‏ ۳۷ بازه رودخانه۱۹ ۰ 
۳ بازه رودخانه۵ YA ANV‏ بازه رودخانه ۲۰ ۹ 
۴ بازه رودخانە۶ ANV‏ ۲۹ جاده ۳۰ 
۱۵ بازه رودخانه۷ ۸ ۳۰ آبراهه AV‏ 
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حدول -Y‏ مفاد پر شاخص‌های ارزیایی عملکرد HecRas Jiu.‏ در مرحله واسنحی و اعتبارسنحی حوزه قره‌جای 


رامیان 
5 قاخص l | als‏ انحراف نسبی شاخص انحراف نسبی شاخص 
سال وقوع رخداد سیلابی انكف D‏ شاخص توافق (O‏ | کارایی نش ساتکلیف (drel) zils‏ 
wy iS‏ 
(Erel)‏ 
yay‏ مرحله وا «VW #44 ۰24۸ a‏ ۳ 
۱۳۹۹ مرحله اعتبار = ۰۷۸ ۰4۹ AO AYY‏ + 


نقشه خطر سیلاب نیز شامل ترکیب نقشه عمق» سرعت و تداوم سیل برای جریان سیلابی با دوره بازگشت 
صدساله می‌باشند. خروجی مدل HEC- RAS‏ عمدتاً برای ag‏ نقشه طغیان و نقشه خطر سیل استفاده می‌شود که 
دامنه تغییرات آن در پنج کلاس تعریف می‌شود (نعیم و همکاران. ۲۰۲۱؛ تبارک و همکاران, ۲۰۲۱) که درجات خطر 
در نقشه پهنه بندی شده (شکل ۲) و فراوانی آن در جدول ۳ ارائه شده است. نتایج نشان می‌دهد که عمده کلاس 
be‏ سیل مربوط به کلاس اول به مساحت ۲۷/۱۰ هکتار OVV)‏ درصد) از پهنه قرار دارد که امتیاز نرمال شده OF‏ 
صفر است که بر خطر بسیار کم عمده منطقه موردمطالعه دلالت دارند. 


حدول ۳- Sle jo‏ خطر سیل در دوره برگشت صد سال در رودخانه قره‌جای 


دوره بازگشت ۱۰۰ سال 
ردیف درجه خطر 
مساحت (هکتار) درصد فراوانی نسبی 
۱ ۱ بسار کم ۰ ow‏ 
YO y‏ کم ۱۳۳۹ ۳۹۹۹۵ 
۳ ۵ متوسط 0.۷ 11.۳ 
“o £‏ زیاد ۱۹۹ AA‏ 
۵ ۱ خیلی زیاد ۲ ۹۱ 
جمع £0.47 1۰ 


سال دوازدهم بررسی قابلیت مدل دو بعدی HEC-RAS‏ در تهیه .... 1۹۹ 


شکل ۲- نقشه خطر سیل برای سیلاب با دوره بازگشت صدساله رودخانه قره‌جای 


۶- جمع‌بندی 

سیلاب فرایندهای طبیعی ناشی از رخدادهای حدی آب و هوایی هستند که Slab‏ شدن دشت‌های سیلابی را 
منجر شده و خسارات اقتصادی» محبط‌زیستی و تلفات انسانی bE‏ توجهی را به بار می‌آورد. آنچه که مورد انتظار 
است این است که تغییر اقلیم بر فراوانی و پهنه سیلاب از طریق SE‏ آن بر شدت و سرعت چرخه هیدرولوژیکی 
تأثیر می‌گذارد (پوینز و تیم ۲۰۱۹). یکی از نقشه‌های مورد نیاز برای ارزیابی ریسک سیل تعیین خطر سیلاب است 
که شامل احتمال وقوع سیل با بزرگی مشخص در یک دوره معین از زمان و مکان مشخص می‌باشد که با بررسی 
عوامل دیگر شامل در معرض قرارگرفتن که از ارزش اقتصادی زندگی انسان و دارایی‌های متأثر از پهنه سیلاب 
حکایت دارد و آسیب‌پذیری که درجه اجزای مختلف (مردم. ساختمان, زیرساخت فعالیت‌های اقتصادی) است 
می‌توان به ارزیابی ریسک سیل یک منطقه دست یافت (فری" و همکاران, ۲۰۲۰). 

در این مطالعه از بررسی‌های عمق سیل رخداد ۱۳۹۷ و رخداد سیل با دوره بازگشت صد سال می‌توان نتیجه 


گرفت که با افزایش دوره باز کشت عمق سپلاب افزایش می‌یابد (پوینز و تیم ۵۸۹ و عمده کلاس‌های خطر با 
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درجه بالا در مناطق میانی از بستر رودخانه در دوره بازگشت صد سال که عمیق‌ترین مناطق را تشسکیل می‌دهند 
قراردارند و حواشسی رودخنه از میزان خطر کمتری برای مردم برخوردار اتتا از نظر مرگ یا جراحته ناچیز 
است و اقطان می‌رود حسارات وارده به اموال نسبتاً کم باشد (دینه" و همکارانء ۲۰۱۲ همچنین با مقایسه پهنه سل 
رخدادهای سال‌های ۱۳۹۲ ۱۳۹۷ و دوره بازگشت صدساله می‌توان اذعان داشت که با افزایش دوره بازگشت مناطق 
مسکونی واقع در روستای سیدکلاته متأثر از سیلاب می‌شوند (مایوسپنتیل " و همکاران» ۲۰۱۹). به‌طورکلی ظرفیت 
کانال برای عبور دبی با دوره بازگشت صد سال تقریباً کافی است ولی در بخش‌هایی از روستای سیدکلاته سیلاب از 
سمت چپ بازه رودخانه وارد روستا می‌شود. به‌طورکلی هدف از نقشه ریسک سیل و ارزیابی آن درک احتمال وقوع 
سیل با شدت خاص که در یک دوره زمانی طولانی رخ می‌دهد است و اینکه بتوان از انجام فعالیت‌های مدیریت 
ریسک که در پهنه سیل حمایت کرد. 

همچنین بررسی polis‏ شاخص‌های ارزیابی مورد استفاده مدل دلالت بر موفق بودن شبیه‌سازی و روندیابی 
جریان سیلابی دارد و این مدل از دقت خوب و مناسبی برخوردار اسست. مدل دوبعدی HEC-RAS‏ جریان 
E peels‏ قل plats Fle‏ سلا lon‏ مس کد ابا خر وس ما فار goss (0 ea‏ 
مشاهداتی نظیر عمق و برآورد سرعت در نقاط مختلف واسنجی شوند. توانایی مدل برای تخمین‌های قابل‌اعتمادتر و 
روندیابی دقیق سیل بستگی به داده‌های مورد استفاده از قبیل DEM‏ مقادیر زبری و تنظیمات اجزای محدود مش" 
aji‏ 

کمبود سوابق و اطلاعات رخدادهای سیلابی گذشته یکی از چالش‌های پیشرو در شبیه‌سازی جریان سیلابی 
است. استفاده از مدل رقومی ارتفاع با قدرت تفکیک مناسب. اندازه مناسب و AS‏ برای مش محاسباتی و نیز لابه 
ضریب زبری در شبیه‌سازی مدل HEC-RAS‏ دوبعدی به‌عنوان مدل al‏ کمک شایانی می‌کند. 

طی مطالعه حاضر مشخص شد حساسیت مقادیر زبری مانینگ در بستر رودخانه در دامنه مناسب در مقیاس 
رخداد سیل WAV‏ زیادتر می‌باشد. یکی از مهم‌ترین پارامترهای هیدرولیکی برای واسنجی مدل مورد استفاده ضریب 
زبری می‌باشد که بر اساس جداول موجود و نظرات کارشناسی که معیار مشخصی برای بهترین مقدار ضریب زبری 
است. انتخاب می‌شود. مدل دوبعدی HEC-RAS‏ از وابستگی کمتری نسبت به‌اندازه ابعاد سلول محاسباتی برخوردار 
بوده و همین امر باعث می‌شود بتوان به‌صرفه هزینه کمتر در زمان کمتر به شبیه‌سازی مناسب جریان پرداخت. 

ازآنجاکه تهیه نقشه خطر سیل با توجه به اطلاعات مهمی چون عمق آب. سرعت. تداوم وسعت پهنه سیل که در 
اعتیار می‌گذارد و نشان‌دهنده شدت سیل و احتمال تجاوز مرتبط با آن‌ها در یک منطقه می‌باشند (دی‌بلداساری و 
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از حمایت‌های دانشگاه لرستان در انجام رساله این‌جانب» که بخشی از آن در مطالعه حاضر ارائه شده است؛ 


صممانه سپاسگزاری می‌شود. 
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